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We present here a simple and generic model for the lumi- 
nous {Q 1 > 10 erg) and long duration (t 7 ~ 10s) Gamma 
Ray Bursts (GRBs) based on the fundamental fact that the 
General Theory of Relativity (GTR) suggests the existence of 
Ultra Compact Objects (UCOs) having surface gravitational 
red-shift z s < 0.615 even when most stringent constraint is 
imposed on the equation of state. This simple model may ex- 
plain the genesis of an electromagnetic fireball (FB) of energy 
as high as Qfb ~ 5 x 10 53 erg and an initial bulk Lorentz 
Factor as high as r\ ~ 10 3 . 
PACS: 98.70. Rz, 97.60.-s 



It is now clear that a large number of Gamma Ray 
Bursts (GRB) involve emission of 7-rays as large as 
Q 7 ~ 10 52 — 10 53 erg under condition of isotropy Q|. Af- 
terglow observations of GRB970228, 970508 and 980703 
show that they indeed have quasi spherical morphology 
H|). In fact, if GRB9901023 were also isotropic, one 
would infer a value of Q 1 w 3.4 x 10 54 erg ||]! However, 
in this paper we shall not consider the unique case of 
GRB9901023, which might be anisotropic jj]], an d focus 
attention on the (other) most luminous events recorded 
so far. We explain them below as events related to the 
formation of UCOs whose existence is suggested by GTR 
irrespective of the details of the EOS of the collapsing 
matter. Since this is a spherical model, unlike the irreg- 
ular non-spherical models, the liberated energy will be 
in the form of photons and neutrinos alone, and, not in 
gravitational radiation. 

As was first shown by Schwarzschild in 1916 ||, GTR 
yields an absolute upper limit on the value of the surface 
gravitational redshift of a static relativistic spherical star 
irrespective of the details of the Equation of State (EOS): 
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Here the subscript "s" refers to the respective "surface" 
values, R is the invariant circumference radius, c is the 
speed of light, and M is the gravitational mass inclosed 
within R = R 
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where p is the total mass-energy density, dV — 4irR 2 dR 
is coordinate volume element, and the symbol dM is self- 
explanatory. Schwarzschild obtained this limit for homo- 
geneous stars by demanding that the central pressure of 



the star does not blow up. This result is actually valid 
even for non-homogeneous stars (see pp.333 of Weinberg) 
H and is obtained when the EOS is allowed to have a 
causality violating sound speed c s = (dp/ dp) 1 / 2 > c. 
When the EOS is constrained to obey causality, it fol- 
lows from Eq.(9.5.19) of Shapiro & Teukolsky (pp.261) 
0, that one would have a tighter limit on z c = 1.22. If 
one constrains the EOS further so as to have c s < c/\/3, 
it follows H that one has an even tighter bound on 
z c = 0.615. To present a realistic model, in the following 
we shall work with the tightest GTR and EOS bound on 
z s = z c = 0.615. It may be also noted that this limit on 
z s is independent on the precise value of M itself. Thus, 
this limit may be applicable to both stellar mass com- 
pact objects like Neutron Stars (NSs) or even supermas- 
sive stars, and, hence, it was debated after the discovery 
of quasars, whether their redshifts were of gravitational 
origin Note that, the presumed canonical NS has a 
value of M ~ 1M Q and R ~ 10km with 2GM/Rc 2 ~ 0.26 
and z s ~ 0.16. However, actually, many existing EOSs 
easily allow a value of M = (2 — 3)Mq and R » 7Km 
(for degenerate matter, there is inverse relationship be- 
tween M and R) [Q, and this may result in a value of 
2GM/Rc 2 « 0.63 or z s w 1.0. Thus, we can very well 
have an UCO, in lieu of a so-called NS. And our UCO is 
nothing but a NS having a compactness, though, higher 
than the canonical (i.e, assumed) value, very much al- 
lowed not only GTR but also by all existing EOSs. So, 
all that GTR tells here is that for a non-trivial (p ^ 0) 
EOS, gravitational collapse process may end with static 
objects having z < z c . Beyond, z > z c , there would be 
no stable static configuration. However, light can escape 
from the collapsing body and the collapse is, in principle, 
reversible, until one crosses a deadline of 2GM/Rc 2 = 1 
or z — 00, when the collapse becomes irreversible. How- 
ever, because of inaccuracies, the present day numerical 
computations have a blurred vision about z c , and further, 
erroneously, they conclude that the collapse becomes ir- 
reversible immediately after z — z c « 0.16. 

Recall that the 10 GTR collapse equations form a set 
of highly complicated coupled non-linear partial differ- 
ential equations. Even, numerically, they can be solved 
to obtain a unqiue result only for the homogeneous dust 
(p = 0). Other cases might also be solved upto a certain 
extent but there could be hundreds of solutions depend- 
ing on a number of explicit or implicit assumptions one 
makes, like self-similarity, adiabaticity, polytropic EOS, 
the variation of the polytropic index, nature of opacities, 
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radiation transport properties etc, etc. Equally impor- 
tant is the question of the initial conditions one explic- 
itly or implicity assumes. And then depending on the 
expectations, one might get the desired result. The gen- 
esis of a high z s object would be marked by emission of 
energy flux Q ~ Mc 2 and then it becomes practically 
impossible to handle the most complex coupled energy 
transport problem in a precise manner. By definition, 
in such cases, one requires to work in the strong grav- 
ity limit where most of the inherent assumptions break 
down. For stellar mass objects, although, the high den- 
sity cold EOS is known with relatively more certainty, 
our knowledge about finite temperature EOS of nuclear 
matter at arbitrary high T is, at present, at is infancy. 
Also note that for numerical computations, for a total 
accumulated uncertainty of few percentage, arising from 
either present theoretical inputs and intrinsic simplifi- 
cations, a potential result like 2GM/Rc 2 rs 8/9 (corre- 
sponding to z s w 2, with finite gravitational acceleration) 
may precipitate to a "2GM/ Rc 2 w 1" (z « oo with in- 
finite gravitational acceleration) signalling the apparent 
formation of an early "event horizon" or a "trapped sur- 
face" . Even if we consider the infinitely simpler problem 
of collapse of an inhomogeneous dust, there could be var- 
ied numerical results, and, in particular, there is a raging 
debate whether such collapse gives rise to a Black Hole or 
a "naked singularity" l|. Such gross uncertainties may, 
at present, obfuscate the signals of formation of more 
compact NSs. In any case, as dicussed before, both the 
nuclear EOSs and GTR actually suggest the existence of 
more compact NSs. 

The self-gravitational energy of a static relativistic star 
is given by [g 
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Then recalling the definition of z from Eq.(l), we may 
write 

E g = - J z{R)c 2 dM w -az s Mc 2 ~ -z s Mc 2 (4) 

where a ~ 1 is a model dependent parameter. The bind- 
ing energy, i.e., the energy liberated in the formation of 
the eventually cold stellar mass compact object, is given 
by virial theorem to be Eg ~ (1/2) I E g |. Most of this 
binding energy is expected to be radiated in the form of 
v — v during the final stages of formation of the UCO: 
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So, given the most restricted limit z c — 0.615 the max- 
imum value of Q v w 0.6M Q c 2 M 2 « 1.2 x 10 54 Af 2 erg 
where of M = M 2 2M Q . This is in agreement with our 
similar previous crude estimate JtJ . The value of Q v mea- 
sured near the compact object will be higher by a fac- 
tor (1 + z s ): Q' v = z s (l + z s )M/2 (now we set c = 1). 



For the NS-formation case, the neutrinos diffuse out of 
the hot core in a time t v < 10s and we may expect a 
somewhat longer time scale for the diffusion of neutri- 
nos from the nascent hot UCO. However, here note that, 
the rather long value of t v < 10 s occurs because of co- 
herent scattering of neutrinos by the heavy (Fe) nuclei 
. If the Fe-nucleons are already partially dissociated by 
an immediately preceding heating, the rise in the value 
of t v for the UCO formation need not be much larger. 
And the locally measured duration of the burst would 
be t' u = (1 + z s )~ 1 t u . Therefore, the mean (local) v — v 
luminosity will be 
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where t v = <iolOs. It may be noted that this value 
of L' v is well below the corresponding j/-Eddington lu- 
minosity [QJ|]. The luminosity in each flavour will be 
L'i = {\/2>)L' v . By assuming the radius of the neutri- 
nosphere to be i?„ « R, the value of effective local neu- 
trino temperature T' (assumed to be same for all the 
flavors), is obtained from the condition 



L' = ^^R 2 aT' A 
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where a is the Stephan-Boltzman constant. Therefore, 
we have, 



T' = 
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where R = RqIO . For a Fermi-Dirac distribution, the 
mean (local) energy of the neutrinos is E' v w 3.15T' ~ 48 
MeV (for z s = .6). The various neutrinos will collide with 
their respective antiparticles to produce electromagnetic 



pairs by the v - 



- e process. The rate of energy 



generation by pair production per unit volume per unit 
time, at a distance r from the center of the star, is given 

by %■ 
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Here, L[(r) 



is the ^-luminosity of a given flavour 



above the ^-sphere, G% = 5.29 x 10~ 44 cm 2 MeV~ 2 
is the universal Fermi weak coupling constant squared, 
K V i = 2.34 for electron neutrinos and has a value of 
0.503 for muon and tau neutrinos. Here the geometrical 
factor ip(r) is 

<p{r) - (1 - xf(x 2 + Ax + 5); x = [1 - (Rjr) 2 ] 1 / 2 

(10) 
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Now, considering all the 3 flavours, a simple numerical 
integration yields the local value of pair luminosity pro- 
duced above the neutrinosphere : 
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This estimate is obtained by assuming rectilinear propa- 
gation of neutrinos near the UCO. Actually, in the strong 
gravitational field near the UCO surface the neutrino or- 
bits will be curved with significant higher effective inter- 
action cross-section. Since, most of the interactions take 
place near the ^-sphere, for a modest range of z s , we may 
tentatively try to incoroprate this nonlinear effect by in- 
serting a (l + z s ) 2 factor in the above expression. On the 
other hand, the value of this electromagnetic luminosity 
measured by a distant observer will be smaller by a fac- 
tor of (l + z s ) -2 , so that eventually, L± = L' ± of Eq.(ll). 
And the total energy of the electromagnetic FB at oo is 
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Thus, the efficiency for conversion of Q v into Qfb is 
Qfb 
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In particular, for z s = 0.615, M 2 = 1, R% = 1 and t w = l, 
we obtain a large e± ~ 15.5%, and it may be reminded 
here that the value of e± should saturate to a limiting 
value of ~ 40%, corresponding to a local statistical equi- 
librium between the 3 flavours of v, v and e + ,e~. This 
highest value of efficiency may be attained, for instance, 
for R — 7km and M — 2.5M Q . Correspondingly, we ob- 
tain a highest value of Qfb ~ 4.8 x 10 53 M2 erg in this 
model. And thus we may explain the energy budget of 
GRB971214, Q 1 w 3 x 10 53 erg 0] without overstretching 
any theory or making any unusual assumption or invok- 
ing any unconfirmed exotic physics (like "strange stars'). 

Now, we shall address the question of baryonic pol- 
lution: 77 = Qfb /AM > 10 2 . In general all mod- 
els involving collision and full/partial disruption of com- 
pact object(s) will spew out thick and massive debris 
(few M Q > M* > O.1M ). Part of this debris is likely 
to settle into a torus and an uncertain small fraction 
(AM) may hang around the system and get accreted 
on a long time scale or may even be unbounded. It is 
practically, impossible to simulate the latter fraction dy- 
namically even in a Newtonian theory. On the other 
hand, spherical implosion models are completely free from 
the presence of such unaccountable and intractable thick 
collisional debris. However, in a normal SN event (as- 
sumed to be basically spherical implosion), the ejection 



of baryonic mass ~ 0.1M Q occurs probably because of 
shock mediated hydrodynamic process. Since, by defi- 
nition, the system is gravitationally bound, any normal 
hydrodynamic attempt of mass ejection can not be much 
successful in a spherical model. But the shock generates 
additional entropy and heat in its vicinity and might be 
able to effect the mass ejection. Yet, the shock is con- 
stantly depleted of energy and gets stalled because of 
^-losses, and disintegration of heavy nuclei ||. Proba- 
bly, the shock might be rejuvenated by the "shock re- 
heating mechanism" B. The energy transfer between 
neutrinos and matter behind the shock is mediated pri- 
marily by the charged current reactions v e + n — > p + e~ 



and v P 
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n + e + . When these reactions proceed to 



the right, the matter heats up, and conversely, the mat- 
ter cools. To have a successful and sufficient net heating 
is a critical phenomenon, and present day (realististic) 
SN codes are unable to find the shock mediated mass- 
ejection (explosion) even in a relatively weak nascent-NS 
gravitational field jl(J. It is not surprising then that 
the same numerical calculations, at present, do not find 
existence of UCO, whose study involves strong gravity, 
finite temperature EOS and complex physics. Probably, 
only for a narrow range of initial conditions and modestly 
deep gravitational potential well this mechanism of shock 
ejection is successful. Thus the real issue is not how to 
explain the non-ejection of mass by direct hydrodynamic 
processes by defying the extremely deep relativistic po- 
tential well. On the contrary, the meaningful question 
is how, for a weak Newtonian potential well, for certain 
range of initial conditions, there could be successful hy- 
drodynamic mass ejection. Note that an UCO with a 
modest value of z s ~ 0.5 has a potential well which is 
~ 300% deeper than the one associated with a canonical 
NS, z s ~ 0.16. Again the basic reason that a critical 
phenomenon like shock heated mass ejection might be 
successful for the SN case is that as one moves from a 
relativistic potential well (high z s ) to a Newtonian well 
(z s < 0.2) is that while the local temperature due to v- 
hating may decrease slowly T" ~ z° 25 and the ^-matter 
interaction cross-section a v<m ~ T' 2 ~ the depth of 
the potential well drops rapidly cx z s . Even then, it is 
far from clear how the hydrodynamic mass ejection can 
really occur. In fact there are ideas that departure from 
spherical symmetry induced by rotation and magnetic 
field might be important in effecting the SN mass-ejection 

I- 

On the other hand, there is a genuine possibility, that 
all models of cosmological GRBs, irrespective of whether 
they explicitly invoke the v + D — > e + + e~ process or not, 
should involve strong direct electromagnetic or v— heated 
mass loss. Even if an unusual pulsar is assumed to emit 
~ 10 52 erg/s rather than 10 38-40 erg/s, the superstrong 
return current impinging back on the pulsar may drive a 
catastrophic wind, a possibility overlooked so far. On the 
other hand, for the thin outermost layers of the object 
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(UCO or an hot accretion torus) emitting the neutrinos, 
well above the j'-sphere, the v- flux S u may induce a 
super-Eddington photon flux S P h • 

Also, though, for a torus with uncertain dynamically 
changing geometry, it is difficult to make any semi- 
analytical or numerical estimate of such a process, in gen- 
eral this effect is expected to be much more pronounced 
because its gravitational self-binding ( z s ) is much weaker 
than that for a spherical UCO surface. And even if a 
steady state model calculation yields a high value of n, 
the eventual value of r\ might be very low if the jet is 
intercepted by this debris. In fact, (by ignoring such 
unmanagable real life uncertainties and difficulties), de- 
tailed Newtonian and crude post Newtonian calculations 
for the NS-NS collision case have been presented by sev- 
eral authors and the conclusion is that, it is difficult 
to understand a value of r\ higher than few. 

For non-spherical configurations, it is difficult to ensure 
that a small fraction of the accreted matter itself is not 
contaminating the FB. And the estimate of AM may be 
made with much larger confidence only for a spherical 
model, where by definition, the entire matter, in general, 
is moving inwardly. Here, a certain fraction of the matter 
lying above the neutrinosphere may be ejected out by the 
^-heating and it may be possible to crudely estimate the 
baryonic mass lying above the i/-sphere independent of 
the details of the problem. The mean cross section for 
i/ e -matter interaction is approximately given as H 
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so that, given our range of E' v ~ 40 — 50 MeV, the value 
of ov.m < 10~ 40 cm 2 . Since the ^-optical thickness of 
the layer above the i/-sphere is ~ 2/3, the surface den- 
sity of this layer S ~ 2m p /(3<7 1 , im ) ~ 10 16 g/cm 2 , where 
nip is the proton rest mass. Therefore, the mass of the 
matter above the ^-sphere is AM ~ AttR^S ~ 10 29 g ~ 
10~ MqRq. Probably, the most detailed work on this 
problem of j/-driven mass ejection from a hot nascent NS 
is due to Duncan, Shapiro & Wasserman jy]; and the 
Table 5 of it shows that for R w 10 6 cm, M = 2M Q , 
we have AM w 1O~ 4 M , if T' = 20MeV. On the other 
hand, for T' = 30MeV, one has, AM « 7 x 10~ 4 M o . 
This estimates were made in the framework of Newto- 
nian gravity, and a GTR calculation, if possible, would 
certainly yield, lower values of AM. Even considering 
these Newtonian values of AM, we find that the value of 
77 could be easily lie between 10 3 > 77 > 10 2 . Now the oc- 
currence of luminous and long GRBs can be understood 
by using the existing ideas Jl3[ . 

Previously, accretion induced collapse of a White 
Dwarf to a NS was suggested as a model for GRBs || . 
The difficulty of this model was that (i) for a canonical 
NS with M = 1M and z s w 0.16, Eq.(12) yields a very 
low value of Qfb ~ 4 x 10 50 erg, (ii) further, the value 
of r) is seen to be too low, and (iii) in the correspond- 



ing weak gravitational well, it can not be ensured that 
supernova shock is not launched. 

And it is probable that less luminous and short GRBs 
may occur by various other scenarios too including the 
so-called "collapsar" or "hypernova" type models JI4J . 
which, may explain the origin of a value of Q 1 ~ \§ — 
10 51 erg for a duration of t~ t < Is (though in extreme 
cases they go beyond this range). 
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